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1 Introduction 



Recent astrophysical data indicate that our universe is currently in a phase of accelerated expansion. 
The physical origin of this acceleration is not completely understood and the related issue is commonly 
called the dark energy problem. 

Several possible explanations have been proposed in the literature. One of them is based on 
the use of modified gravitational models, the simplest one consisting in the inclusion of a small and 
positive cosmological constant. Such model works quite well, according to the most recent data, but 
nevertheless it has some drawbacks (see, for example [UEKSlll] and reference therein). 

Roughly speaking, the idea underlying modified gravity models is that the Einstein-Hilbert action 
gives only an approximate low energy contribution to gravitation and additional terms depending 
on the quadratic curvature invariants should necessarily be included. This idea is quite old since it 
was already contained in the seminal paper [S], where quadratic terms in the curvature, justified by 
quantum effects, were added to the Einstein-Hilbert Lagrangian (for a review, see [6]). The important 
recent finding is that the inclusion of suitable higher order contributions may realize not only the 
actual accelerated expansion, but also the early time inflation epoch j^. Within this context, the 
de Sitter (dS) space-time plays a fundamental role, being able to provide an acceleration at different 
stages in the cosmological set up. 

In previous papers |8i [9l [TOl [11] , f{R) gravity models and a non local Gauss-Bonnet gravity model 
at one-loop level in a de Sitter background have been investigated. A similar program for the case of 
pure Einstein gravity was initiated in Refs. [121 [131 IH] (see also [151 [16]). Furthermore, such approach 
also suggests a possible way of understanding the cosmological constant issue |14j . Hence, the study 
of one-loop generalized modified gravity is a natural step to be undertaken for the completion of such 
program, with the aim to better understand the role and the origin of quadratic corrections in the 
curvature. An alternative approach, which is in some sense alternative, has been proposed by Rueter 
and collaborators [17J, see also the review paper [181, and [19], in which quantum gravity effects in 
astrophysics and cosmology are presented. 

In the present paper, we will investigate in some detail a model described by a Lagrangian density 
which depends on geometric quadratic invariants. The quantization of quadratic models of gravity 
has been discussed in many papers, and in particular studied in detail on fiat space in the seminal 
paper |20j . A preliminary discussion of a quadratic model based on one-loop on-shell results has been 
presented in [2T] . 

Here we start with the classical Euclidean gravitational action 

lE[g] = -j (fx^F{R,P,Q) = -j (fxy/g[M^R-2a + b{R^-?,P) + pG], (1.1) 

where h and j3 are dimensionless parameters, M"^ a mass-squared parameter playing the role of gravi- 
tational coupling constant, and a a "cosmological constant" dimensional term. By G we indicate the 
Gauss-Bonnet topological invariant which, in 4-dimensions, does not contribute to the classical field 
equations. For this reason, the action in (jl.ip in classically equivalent to the so called Einstein- Weyl 
gravity, since the quadratic Weyl invariant W and the Gauss-Bonnet invariant G are related by 

G-W = \{R^ -W), G = R^-AP + Q, W = \r'^-2P + Q, (1.2) 

P and Q being 

P = R'^Rij, Q = R'^'^'Rijrs, i,i,r,s = 0,1,2,3. (1.3) 

As we already said above, at classical level the Gauss-Bonnet term does not play any role and could 
be dropped off but, as we shall see in the following, it do will play an important role at quantum level. 

The action in the form (jl.ip is quite useful in order to discuss the so called "critical gravity", 
which corresponds to a particular choice of the b parameter. An extensive study of 4-dimensional 
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critical gravity, in the presence of a negative cosmological constant, has been recently presented in 
Refs. \22 \ \23\ [2 ^ I25j. where additional relevant references can be found. 

It should be noted that one-loop Euclidean quantum gravity in a de Sitter background — as was 
fully exploited in [T3| for the case of Einstein's gravity in the presence of a cosmological term — presents 
some peculiar aspects within the background field method. First, working with the Euclidean version 
5(4), one is dealing with a geometric background associated with a compact manifold without a 
boundary. This means that the volume is finite and can be expressed as a function of the constant 
Ricci curvature, which may be chosen as background field. A second important remark is that, in 
order to discuss the one-loop renormalizability of the model, as well as the related renormalization 
group equations, one is forced to work with the off-shell one-loop effective action. As a consequence, 
the Landau gauge appears to be the most convenient one. Besides, the usual effective action calculated 
in this gauge coincides with the Vilkovisky-De Witt effective action (see, for example, [[Bj). 

Such approach should be compared with the more traditional one, nicely reviewed in [26], where 
the Sakharov induced gravity approach [27] and its modern variants [25] have been discussed too. 
Conformal gravity has been discussed in [29]. The ghost absence issue for a very general gravita- 
tional quadratic model on Minkowski space-time has been recently investigated in [30]. Furthermore, 
alternative approach is presented in |31] . 

Regarding to the choice of regularization, since we are dealing with non flat space-time, it is almost 
mandatory (or at least very convenient) to make use of a variant of the generalized zeta-function 
regularization [32l [33] (see also [331 ESI |36] ) , and the associated heat-kernel techniques [371 ESI [39] . In 
this way, one may evaluate the one-loop effective action and then study the possibility of stabilization 
of the de Sitter background by quantum effects. 

The paper is organized as follows. Section II contains the evaluation of the quantum fluctuation 
operators relevant for the one-loop calculations to be carried out. In Section III, the off-shell one loop 
partition function is presented and the corresponding one loop renormalization is discussed. Finally, 
Section IV is devoted to conclusions. 



2 Quantum field fluctuations around maximally symmetric instan- 
tons 

In this Section we will discuss the one- loop quantization of the model in (jl.ip on a maximally symmetric 
space (see, for instance [B]). To start with, we consider the Euclidean gravitational action in 
and, for convenience, we separate linear and quadratic terms 

F{R,P,Q) = f{R) + b{R^ -3P) + /3G, f{R) = M^R-2a. (2.1) 

The model admits a constant Ricci curvature solution Rq. In fact, the general equation for the 
existence of de Sitter solution [lOl HI] 



is trivially satisfied, and reads 



= 0, (2.2) 

R=Ro 



1 Aa 
f{Ro)--Rof'{R^) = =^ ^' = J^- (^-^^ 

We are interested in studying quantum fluctuations around the Euclidean dS instanton S"^ with 
positive constant scalar curvature Rq. This is a maximally symmetric space having covariant conserved 
curvature tensors. Its metric may be written in the form 

dsl = dr\l - Hy) + (Y^^ + r'dSl , (2.4) 
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dS'^ being the metric of the two-dimensional sphere S'^ and Hq the Hubble constant. The finite volume 
is given by 



3847r^ 
Rl 



Rq — 12Hq , 



Gn 



while the Riemann and Ricci tensors are 



R^ 
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ijrs 
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Now let us consider small fluctuations around the maximally symmetric instanton. For the sake of 
completeness, we consider the general action discussed in [2T], but linear in P,Q. Then, we shall 
restrict to the action (|l.ip . at the end of the computation. For simplicity, we also put = 1. When 
necessary the right units will be easily recovered by dimensional analysis. 
We set 



9ij ^ 9ij ~l~ hij , 



where from now on gij = gf^^ is the metric of the maximally symmetric space and, as usual, indices 
are lowered and raised by means of such metric. 
Up to second order in /i^, one has 



_ + h'^h{ + 0{h-^), h = g'^ hij , 



(2.7) 



V9 — > V9 



and 



R 



Rq 
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^ h + V,Vjh'^ - Ah 



+^ h^^hjk - 7 VihV'h - - VkhijVh'^ + VihlVjh^' 
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V^hikV'h^' 



(2.9) 



where Vfc represents the covariant derivative in the unperturbed metric gij . More complicated expres- 
sions are obtained for the other invariants P, Q, but for our aims it is not necessary to write them 
explicitly. 

By performing a Taylor expansion of the Lagrangian around de Sitter metric, up to second order 



in hij, we get 



F{Ro,Po,Qo) + ^ + C2 



(2.10) 



where £2 represents the second-order contribution and X vanishes when the de Sitter existence con- 
dition ([221) is satisfied. For our particular model, X = [f{Ro) - {l/2)Rof'{Ro)]/M'^. 

It is convenient to carry out the standard expansion of the tensor field hij in irreducible components 
[14| . namely 



hij = hij + Viij + Vjii + ViVja + - gij{h - Aa) , 



(2.11) 



where a is the scalar component, while (,i and hij are the vector and tensor components, with the 
following properties 



v^r = 



Vih'^ = 



h] = 0. 



(2.12) 



In terms of the irreducible components of the hij field, the Lagrangian density, disregarding total 
derivatives, becomes 



Co 



C-hh + 2 Cha + + Cy + Ct ; 



(2.13) 
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where Chh-, ^^ha-, ^^tra represent the scalar contribution (a 2 x 2 matrix), while Cy and Ct represent 
the vector and tensor contributions, respectively. One has 

1 



C 



hh 



h 



32 



FrrK 



|-F^A + ^FpRoA + ^FqRoA 
62 lb lb 



+^FrrRoA + ^FpA^ + ^FqA^ + ^FrrA"" 
lb lb lb 62 



h, 



(2.14) 



-—FrrR'oA + —FrRoA - -FpRoA' 
lb lb 

-^FqRoA^ - ^FrrRoA^ + ^FrA^ 



-FpA' 



^FnA^ - -FrrA' 
g y 16 



(2.15) 
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16 



+i^FpA'' + ^FqA" + ^FrrA' 
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(2.17) 



— FpRi + —FnRn FrRo --X + -FrA 

72 ° 36^° 24 4 4 

+^Fpi?oA - ^FqR^A + ^FpA + FqA^ . 



(2.18) 



where A = g'^^Vi^j is the Laplace-Beltrami operator in the unperturbed metric gij, which is a 
solution of the field equations, but only if X = 0. We have written the above expansions around a 
maximally symmetric space, which in principle would not be a solution. This means, in other words, 
that the function f{R) can be arbitrary. In the latter expression Fr,Frr represent the first and 
second derivatives of F{R,P,Q) with respect to R evaluated on de Sitter metric Qij. And similarly 
for Fp,Fq. 

As is well known, invariance under diffeomorphisms renders the operator in the {h, a) sector not 
invertible. One needs a gauge fixing term and a corresponding ghost compensating term. Here we 
choose the harmonic gauge, that is 



Xj 



-Vihi-^Vjh = 0, 



and the gauge fixing term 
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The corresponding ghost Lagrangian reads [B] 
6x'' 



gh 



6ei 



■a 



(2.19) 



(2.20) 



(2.21) 
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where and Bj^ are the ghost and anti-ghost vector fields respectively, while 5 is the variation of 
the gauge condition due to an infinitesimal gauge transformation of the field. In this case, it reads 



S hi 



Hj — + Vj-ej =^ 
Neglecting total derivatives, one has 



Qij A + Rij . 



C 



gh 



A 



Cfc 



In irreducible components one finally obtains 



-gf 
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Rq 
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where ghost irreducible components are defined by 



Ac 



Bk + Vkb, 



VkB'' 



0. 



(2.22) 
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3 OfF-shell one-loop effective action 

In order to compute the one-loop contributions to the effective action one has to consider the path 
integral for the bilinear part, £ = £2 + + C,gh-, of the total Lagrangian and take into account the 
Jacobian due to the change of variables with respect to the original ones. In this way, one gets [14^ [B] 

(detGy)-^/' 1 D[hij]D[Ck\D[B''] exp j d^x^C^ 

(det Gij)~^'^ det Jf ^ det j]'^ 
X / D[h]D[hij]D[e]D[a]D[Ck]D[B'']D[c]D[h] exp 



d x^C 



(3.1) 



where Ji and J2 are the Jacobians coming from the change of variables in the ghost and tensor sectors, 
respectively They read 



Ji = A, 



J2 



A 



Ao + ^l A 



(3.2) 



and the determinant of the operator Gij is trivial in this case. Here and in the following Aq, A 1, A2, 
represent the Laplacian acting on scalars, vectors and tensors, respectively. 

Due to the presence of curvature, the Euclidean gravitational action is not bounded from below, 
because arbitrary negative contributions can be induced on R by conformal rescaling of the metric. 
For this reason we have also used the Hawking prescription of integrating over imaginary scalar fields. 
Furthermore, the problem of the presence of additional zero modes introduced by the decomposition 
()2.1ip can be treated by making use of the method presented in Ref. [14J . 

Now, for the action p.ip a straightforward computation leads to the off-shell one-loop contribution 
to the "partition function". In the Landau gauge, /9 = 1,7 — >• 00, with X = Ro/2 — 2a/M^, we get 



of f— shell 



(1) 
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(3.3) 
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of f— shell 



2^- logdet 



^logdet (-1 
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2a 
M2 
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-i logdet (i, [-A 



y+1 ) + ^ logdet f-^ 



-A. 



-An-^ 



where 




2M^ 1 I \ 

-3Ro — ± -^96ba + 4M4 - 20bM^Ro + PR^ 1 . 



(3.4) 



(3.5) 



As usual, an arbitrary renormalization parameter 1/fJ? has been introduced for dimensional reasons. 
As expected, the parameter /3 does not appear in the latter expression, since the Gauss-Bonnet in- 
variant does not give contributions to the field equations, but it gives a constant contribution to the 
classical action, which will actually play an important role in the renormalization procedure. 



3.1 On-shell one-loop effective action 

As is well known, the on-shell effective action does not have to depend on the gauge and, in fact, 
setting X = 0, that is M'^Rq — 4a = 0, we get 



on— shell 



+ 



967r' 
1 



2M^ 
Ro 



1 



+ b] + 64tt^P - - log det 



1 



log det 



-A2 + 



Rq 
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+ ^ log det ( 
2 \ /i^ 



(3.6) 



The above expression is only formal, and one needs regularization. For the moment, let us imagine to 
be dealing with the finite part of such an effective action. 

We observe that there exists a "critical" value for b for which all spin excitations become "massless" . 
In fact, choosing 



-'crit 



2M2 

Rn 



the effective action simplifies to 



crit 



-A 



Rq 

T 



+ log det 



-A2 + 



Rq 



(3.7) 



(3.8) 



In contrast to the AdS case, in Euclidean dS space Rq > 0, and so bcrit < 0. 

The stability of the dS solution can be investigated by looking at the spectra of the Laplace-type 
operators and it then follows that all eingenvalues are non-negative as in general relativity, with the 
possible presence of a zero mode [lUIH]. As a consequence, dS background space is stable, in agreement 
with the classical analysis presented in [301 E] ■ 

But what about the stability of the critical values with respect to renormalization? According 
to the background field method one should work at the off-shell level. Nevertheless, one may try 
an on-shell, one-loop renormalization, by observing that /3 might become a "bare" constant, and its 
redefinition may contain all counterterms necessary in order to cancel the on-shell one-loop divergences 
coming from the functional determinants. In general, using a variant of the zeta function regularization 
procedure [21], at one-loop level one has 



+ C{0\Li)logfi^ + C\0\Li 



(3.9) 
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where the summation is over ah Laplace-type operators appearing in the one-loop contribution to the 
action. 

For the critical gravity in (|3.8|) . we have Li = — A i — i?o/4 and L2 = — A 2 + ^0/6, thus 



r(/^,e) = 64^2(^ + /3(^)) + ^[C(0|Li)-2C(0|L2)] 

+ log 11 [C(0|Li) - 2C(0|L2)] + \ C'(0|Li) - C'(0|L2) , 



(3.10) 



where /3o is the bare coupling constant and /3(/u) the running one. Making a suitable choice for /3o, 
one has the renormalized on-shell effective critical action 



T,rM = 647r2/3(^) +log/i [C(0|Li) - 2C(0|L2)] + ^C'(0|Li) - C'(0|i^2) • 



(3.11) 



The eigenvalues of the Laplace type operators on SO {A) are well known and in this way it is possible 
to compute the zeta- functions appearing in the expression above explicitly. In particular, C(0|Li) = 
—191/30 and C(0|L2) = 89/9, while C'(0|ii) and C'(0|-^^2) are computable expressions independent on 
/i. The usual imposition 



dr(/i) „ 
— = 0, 

gives rise to the renormalization group equation for critical gravity, in the form 
dm 



2C(0|L2) -C(0|ii) ~ 26 > 0. 



(3.12) 



(3.13) 



This is the only running coupling constant and, thus, on-shell critical gravity seems to be stable at 
the one-loop level. But this is not really conclusive since, strictly, the issue of criticality depends on 
the on-shell expression. 



3.2 Off-shell one-loop renormalization 

As far as the off-shell one-loop renormalization is concerned, the situation is completely different with 
respect to the previous one and apparently there is no room for the notion of criticality. 

Again, the starting point is the equation in (j3.9|) . but now the classical action contains all the bare 
quantities, which generate the counterterms for absorbing the one-loop divergencies. It reads 



3847r^ 



3847r' 



M2(//,e) 2a(/x,e) b{n,e) /3{p,,e) 



Ro 



Rl 



+ 



+ 



Ro 



Rl 



+ 



+ 



6 



+ 



(3.14) 



where we have separated the finite and divergent parts of the coupling constants by means of suitable 
finite quantities Ai, Bi, Ci independent of Rq. On the other hand, a direct computation shows that 



Ro 



Rl 



(3.15) 



where Li are all Laplace-type operators in (j3.4p and A{fi) , B (n) , C (fi) are finite functions depending 
on the renormalized running coupling constants M^(/i), Q(/i), 6(/i). They read 



^, , 20M4(u) main) ASa^iu) 



(/i) 



b{^i) M4(/,) ' 



1763 
90 



(3.16) 
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The model is one-loop renormalizable since all one-loop divergences can actually be absorbed by an 
appropriate choice of Ai, Bi,Ci. The finite, renormalized one- loop effective action reads 



M^(/i) 2a(/x) , bin) 

Rq 



Rl 



+ 



+ log^ 



Rq 



Rl 



+ Z, 



(3.17) 



where we have dropped the parameter because here it does not play any role, and we have set 
Z = ^'(0|Li). This is the finite part of the functional determinant which does not depend on ji 

and in principle can be explicitly evaluated. 

As above, the one-loop renormalization group equations can be obtained by means of (|3.12p . To 
this aim it is convenient to introduce the dimensionless variable p = log^, being a reference low 
energy scale. /,From ()3.12p we obtain the three differential equations 

db 



dp ~ '^b , 
dM^ _ 

dp 
da 
dp 



] 

' 1927 



1 / 



Cb 



l2^ + 20ab 



1763 
86407r2 



(3.18) 



where all coupling constants are functions of p. Solving the system of differential equations above, we 
finally get 



b{p) = Cbp + co, 

M\p) = ciip + coy'-P^ [(/o + co)i°P2+C2]'/' 

Pi 



a{p) = 487r 



2 M*ip) 



(p+co) 



„ J 2|32 

+ l+C2{p+Co)-"'f2 



Pi 



[ P2 



863 
17630 ' 

V474769 
17630 



(3.19) 



The integration constants co,ci,C2 depend on the initial conditions and we assume all of them to be 
non negative. Moreover, to simplify the discussion, from now on we shall take C2 = 0. With this 
assumption we get 



M\p) = Clip + cor , 



a{p) = AS-K^p 



oM\p) 



P + Cq 

and the one- loop, running, gravitational coupling constant reads 
1 1 



p = Pi + P2 ^ 0.09 . 



G{p) 



167rM2(/)) 167rci(p + co)P ' 
while the one-loop, running, cosmological constant is 



Hp) 



A'P{p) 



487r pci{p + Co) 



(3.20) 



(3.21) 



(3.22) 



As a result, there is no Landau pole, and at large energy p^ po 01 equivalently p — )• 00, one has that 
both G{p) and A{p) go to zero. This property is the analogue on de Sitter space of the well known 
gravitational asymptotic freedom for quadratic gravity [121 113 IIH 113 SS] ■ 

Furthermore, we may assume general relativity to be valid at low energy, that is p ^ po or 
equivalently /O ~ 0, then 



Gip)\ 



p^O 



Gn 



ci 



167r Gn 



M\p) 



1 



167r Gn 



1 + ^ 

Co 



(3.23) 



Gn = G{0) being the Newton constant. 

To conclude this section, we write down the effective field equation given by 



dr 

dRn 



0. 



(3.24) 
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The solution can be written in the impHcit form 

?3 



Rn 



1 - "^.P 
^ 3847r2M2 



4a Bp ^ dZ 



M2 1927r2M2 3847r2M2 dRo 



(3.25) 



This is a quite comphcated expression in the unknown variable Rq. Of course, at low energy p ~ 0, 
Rq <^ Af2, one gets the classical solution 



^0~T?277TV = 4A(0), (3.26) 



4a(0) 
M|(0) 

but in principle other regimes can be studied 



4 Conclusions 

In this paper, the Einstein gravity plus a quadratic gravitational Weyl term has been investigated by 
computing the corresponding one-loop quantum corrections by means of the background field method. 
As a classical background we have considered the de Sitter one for its potentially very important 
physical applications. The one-loop calculation has been performed in the Euclidean sector, where the 
classical background is the compact manifold S'^. In the calculation, due to the fact that we are working 
with a non flat background, we are forced to make use of (a variant of) zeta-function regularization. 
In the presence of this compact curved manifold, the one-loop effective action, and also the ensuing 
one-loop renormalization group equations, have been computed and carefully investigated. On shell, 
the associated quadratic critical gravity has been discussed. In order to investigate the role of the 
one-loop corrections, we have to consider the off-shell one-loop effective action and so, to get rid 
of the gauge dependence, the Landau gauge has been used. In this way the critical conditions are 
lost, in general. In fact, also in the simplified case we have considered (c2 = 0), the critical ratio 

'^^I'p-j ~ (/O + Co), which is not equal to b{p). This means that the critical gravity condition are not 
stable under the one-loop renormalization flow. As a consequence, one might doubt about its relevance 
at least when the background is the compact manifold S"^. In fact, in the anti de Sitter case (AdS) the 
situation could be completely different, because the Euclidean counterpart of AdS is the non compact 
hyperbolic manifold H^. 

In conclusion, it should also be interesting to repeat this one-loop calculation in alternative ex- 
tended gravity models, for example, in the so-called f{T) gravity models, which depend on geometric 
invariants built up by using torsion field (for details see [17] and the references therein). 
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